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The effectiveness of the perturbative renormalization group approach at fixed space dimension d in the theory 
of critical phenomena is analyzed. Three models are considered: the 0{N) model, the cubic model and the 
antiferromagnetic model defined on the stacked triangular lattice. We consider all models at fixed d = 3 and 
analyze the resummation procedures currently used to compute the critical exponents. We first show that, 
for the 0{N) model, the resummation does not eliminate all non-physical (spurious) fixed points (FPs). Then 
the dependence of spurious as well as of the Wilson-Fisher FPs on the resummation parameters is carefully 
studied. The critical exponents at the Wilson-Fisher FP show a weak dependence on the resummation param- 
eters. On the contrary the exponents at the spurious FP as well as its very existence are strongly dependent 
on these parameters. For the cubic model, a new stable FP is found and its properties depend also strongly 
on the resummation parameters. It appears to be spurious, as expected. As for the frustrated models, there 
are two cases depending on the value of the number of spin components. When N is greater than a critical 
value Nc, the stable FP shows common characteristic with the Wilson-Fisher FP. On the contrary, for N < Nc^ 
the results obtained at the stable FP are similar to those obtained at the spurious FPs of the oIn) and cubic 
models. We conclude from this analysis that the stable FP found for Af < A^c in frustrated models is spurious. 
Since Afc > 3, we conclude that the transitions for XY and Heisenberg frustrated magnets are of first order. 
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1. Introduction 

Renormalization group (RG) methods [l| have led these last forty years to such a deep un- 
derstanding of critical phenomena that it is hard to imagine describing them now without having 
recourse to these methods. In particular, quantitative computations of critical exponents and uni- 
versal amplitude ratios have been achieved this way The prototypical models where these 
methods have been particularly successful are the three-dimensional (3d) 0(7V)-symmetric sys- 
tems modeled by a cj)'^ theory that involves only one coupling constant [3|. Such systems are for 
instance polymer chains in a good solvent (TV = 0) 0|, pure ferromagnets with uniaxial anisotropics 
{N — 1), superfluidity in He'' {N = 2), pure Heisenberg ferromagnets (TV = 3) etc. 

The same approach has been used in the analysis of 3c? anisotropic 0, Q , structurally disordered 
[1, 0, Q ^iid frustrated systems (see also review in [9.]). Formally, these are described by 0^-like 
theories with several coupling constants. There, besides universal exponents and amplitude ratios, 
a subject of interest is the crossover between different universality classes and (universal) marginal 
order parameter dimensions that govern this crossover (lol. 

The prevailing majority of works in this direction aims at getting an accurate description of 
well-established phenomena. As an example, the scalar 0^ theory is used to obtain precise values 
of physical quantities characterizing the critical behavior of models belonging to the 3d Ising 
universality class. The very existence of a second order phase transition in this model is a well 
established fact, observed in experiments and MC simulations and proven by analytic tools. On 
the other hand, there are systems for which the very existence of a second order phase transition 
is not well established and RG can be used to study this point. This is in particular the case 
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when other techniques do not ahow us to conclude. Technically, if RG transformations possess a 
fixed point (FP) which is stable once the temperature has been tuned to T = Tc and whose basin 
of attraction contains the microscopic Hamiltonian, then the transition described by this FP is 
of second order. However, as we show in the following, a FP can be found that has no physical 
meaning and which is only mathematical artifact. In this case, of course, any conclusion based on 
this FP is meaningless. We call it a spurious FP. 

As for XY {N = 2) and Heisenberg [N = 3) frustrated magnets, a stable FP is found within 
the field-theoretical perturbative approach at fixed d [.12- 15J which is usually interpreted as the 
existence of a second order phase transition. However, this is confirmed neither by perturbative s = 
d — A [iB - flst and pscudo-e-expansions fll], nor by the non-perturbative renormalization group 
(NPRG) approach [9, 19-21J. The arguments of [141. where such a FP is found, were questioned 
recently within minimal subtraction scheme (22l. l23l|. In this paper we continue the analysis of 
the non-trivial FPs appearing in the frustrated model along the way started in (23- 24 1. Working 
within the massive renormalization scheme we consider this problem in the general context of FPs 
found within the fixed d approach applied to (/)^-like theories. We use the resummation method 
exploited in [25'| and we study the physical and spurious FPs of the 0{N) models thus showing 
the characteristics of a spurious FP. By comparison, this allows us to classify the FPs obtained in 
the frustrated models. A detailed analysis of 0{N) and frustrated models was performed in d = 2 
elsewhere [2^. Here, only the results for the d = 3 case are reported and they are completed by 
results for the id cubic model. 

The paper is organized as follows. In the next section [5] we present our method of analysis 
of RG functions at fixed d. It is applied to the 0{N) models in c? = 3 in section [31 describing 
stability of 0{N) FPs with respect to variations of the resummation parameters and demonstrating 
differences between physical and spurious FPs. Then the cubic model is investigated in section HI 
where results for the stability exponents at the cubic physical and spurious FPs are presented. In 
section [S] frustrated models with N — 8, N — 2 and = 3 are analyzed. We conclude in sectionjS] 



2. Perturbative RG at fixed d 

Let us start our analysis by introducing the main quantities that are used within the field 
theoretical description of the critical behavior of the 0{N) models. The Hamiltonian reads: 




where is a A^-component vector field and mq > is the bare coupling. 

It is well known that this theory suffers from ultraviolet divergences. Within the field-theoretical 
RG approach their removal is achieved by an appropriate renormalization procedure followed by a 
controlled rearrangement of the perturbative series for RG functions In particular, the change 
of couplings under renormalization is obtained from the /3-function, calculated as a perturbative 
series in the renormalized coupling. The explicit form of this function depends on the renormal- 
ization scheme. However, universal quantities such as critical exponents depend neither on the 
regularization nor on the renormalization scheme at least if the series expansion is not truncated 
at a finite order. A definite scheme must of course be chosen to perform actual calculations. Among 
them the dimensional regularization and the minimal subtraction scheme [26'| as well as fixed di- 
mension renormalization at zero external momenta and non-zero mass (a massive RG scheme) (27| 
are the most used ones. 

Introducing a fiow parameter £ (which can be the renormalized mass) the change of the renor- 
malized coupling constant u under the RG transformations is by definition given by the equation: 

e±uii) ^ (3^{uie)). (2.2) 

A fixed point u* of the differential equation (|2.2p is determined by a root of the following equation: 

f3u{u*)=0- (2.3) 
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The FP is said to be stable if it attracts the RG flow when £ ^ 0. At the stable FP 

(2.4) 



d/3u 
du 



has a positive real part. 

The analysis of the perturbative /3-function can be performed in two complementary ways. The 
first way consists in obtaining u* , the FP values of the coupling, in the form of an expansion in a 
small parameter which is usually e = A — d [28| . In the minimal subtraction scheme, e-expansion 
arises naturally, since e enters the /3-functions only once as the zero-loop contribution. For the 
0(N) models and within the minimal subtraction scheme the /3-function is known at five loops 
(29| . Its two- loop expression is [30|: 

^..(»)-«(—^|^«^ + ■■■)■ <") 

This function has two roots given by a series expansion in e: it* = which corresponds to the 
Gaussian FP and the non-trivial Wilson-Fisher FP: 

u*=e + 0{e^) (2.6) 

which is stable for e > and thus governs the critical behavior of the model for d < A. 

For the massive scheme, d enters all loop integrals and the e-expansion, although possible in 
principle, is in practice performed only at low orders. The function (3u{u) is known in this case at 
six loops [31i] and we give here its two-loop approximation at d = 3 [30|] : 

Instead of the e-expansion, the pseudo-e expansion is widely used in the massive scheme (s^ . It 
consists in replacing the unity in the r.h.s. of equation (12. 7p by the pseudo-e expansion parameter 
T and in looking for a root of as a series expansion in t. Of course, r must be set equal to one 
at the end of the calculation. The non-trivial root u* writes: 

u*^t + 0{t'^). (2.8) 

This technique avoids a situation when errors coming from the solution of an equation for u* 
are included into the series for critical exponents. Therefore final errors accumulate those coming 
from series for /3„(m) and series for critical exponents. Instead, the pseudo-e expansion results in 
a self-consistent collection of contributions for the different steps of calculation [33] . Two common 
features of the e- and pseudo-e-expansions are (i) once the FP is found at one loop, it persists at 
all loop orders (ii) in the limit e 0, all FPs coincide with the Gaussian FP which controls the 



critical behavior of (j)'^ models in d = 4 (see [34ll35| and reference therein) 



The second way to analyze the perturbative /3 functions consists in fixing d and then in numer- 
ically solving equations (|2.3p . This method is usually called in the literature the fixed-d approach 
(27l . [3^ . No real root of can be a priori discarded in this approach contrarily to the e- and 
pseudo-e expansions, where the only FP retained is by definition such that u* e or t [33 ■ As 
a result, the generic situation is that the number of FPs as well as their stability vary with the 
loop-order: at a given order, there can exist several real and stable FPs or none instead of a single 
one. This artifact of the fixed-d approach is already known and was first reported for the massive 
scheme in d = 3 The way to deal with it is also known: the perturbative /3-function has to 
be resummed (see e.g. [s^ and [2|)- Resummation is supposed both to restore the Wilson-Fisher 
FP when it does not exist and to eliminate the spurious ones. The ability of the resummation 
procedures to do so is usually not questioned. However, as we show below on the example of the 
0{N) models, and then for more general models, spurious FPs can exist even after resummations 
have been performed. Some extra criteria are therefore necessary to eliminate them. We describe 
these criteria below. 
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In the minimal subtraction scheme the fixed-c? method was introduced in where e is kept 
fixed to 1 and for the massive scheme calculations were performed in d = 2 and = 3 in (27l . [33 | . 
We note here, that the fixed-c? approach does not necessarily mean that the space dimension d is 
integer: in the minimal subtraction scheme one can perform calculations for fixed non-integer e as 
well as in the massive renormalization one can evaluate loop integrals for non- integer d [40|]. We 
explain in the next subsection our resummation procedure. 

2.1. Resummation 

To extract reliable data from /3-functions one has to resum them. Here, following we use 
the conformal mapping transform. The main steps are the following. For the series 

f{u) = J2anU^ (2.9) 

n 

with factorially growing coefficients an one defines Borel-Leroy image by: 

where r[. . .] is the Gamma-function. This image is supposed to converge, in the complex plane, 
on the disk of radius 1/a where u = — 1/a is the singularity of B{u) nearest to the origin. Then, 
using the integral representation of T[n + b + 1], f{u) is rewritten as: 

oo 

/(«) = y ^ — f dt 6-* ^"+^ (2.11) 



Then, interchanging summation and integration, one can define the Borel transform of / as: 

oo 

/b(u) ^ I dt e"* t'' B{ut). (2.12) 








Figure 1. Conformal mapping of the cut-plane onto a disc. See the text for details. 



In order to compute the integral in (|2.12p on the whole real positive semi-axis one needs an 
analytic continuation of B{t). It may be achieved by several methods. In particular, instead of the 
series of B{t) its Pade-approximants can be used [iilill,!!!!. We prefer to use here a conformal 
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mapping technique. In this method assuming that all the singularities of B{u) lie on the negative 
real axis and that B(u) is analytic in the whole complex plane excluding the cut from — 1/a to — oo, 
one can perform in B the change of variable w(u) = <J=> u(w) — - s'> , as it is shown 

in the figure [TJ This change maps the complex u-plane cut from u = — 1/a to — oo onto the unit 
circle in the ly-plane such that the singularities of B(u) lying on the negative axis now lie on the 
boundary of the circle \w\ = 1. Then, the resulting expression of B{w{u)) has to be re-expanded 
in powers of w{u). Finally, the resumed expression of the series / writes: 

oo 

/r(w) = ^) / H^tT , (2.13) 



where c?„(a, b) are the coefhcients of the re-expansion of B{u{w)) in the powers of w{u). 

It is moreover interesting to generalize the above expression (|2.13p in the following way (4^ 



since this allows to impose the strong coupling behavior of the series: f{u — > oo) ^ u"^^. 

If an infinite number of terms for f{u) were known, expression (j2.14l) would be independent 
of the parameters a, b and a. However, once a truncation of the series is performed /r(u) starts 
acquiring a dependence on these parameters. In principle, a and b are fixed by the large order 
behavior of the coefficients a„ in series (|2.9p [33| . while a is determined by the strong coupling 
behavior of the initial series. However in many cases, only a is known and a and b must be 
considered either as free or variational parameters. In any case, the choice of values of a, a and b 
must be validated a posteriori by checking that a small change of their values does not yield strong 
variations of the quantities under study. 

The above described procedure may be generalized to the case of several variables. For instance, 
when / is a function of two variables u and v, the resummation technique used in [14] treats / as 
a function of u and of the ratio z — v/u: 

f{u,z) = Y,an{z)u^. (2.15) 

n 

Then, keeping z fixed and performing the resummation only in u according to the steps (|2.10p - 
(ETH), one gets: 

Uu, z) = J2 dnia, a{z), 6; z) J dt e'* [^K^)]" (2.16) 
„ J [1 - w{ut; z)] 

with: 



^l + a{z)u-l 

w(u:z)^ — , = (2.17 

^l + a{z)u+l 

where, as above, the coefficients dn{o:, a{z),b, z) in (|2.16p are computed so that the re-expansion 
of the right hand side of (|2.16|) in powers of u coincides with that of (j2.15p . 



2.2. Principles of convergence for results obtained with resummation 

The dependence of the critical exponents upon the parameters a, b and a is an indicator of the 
(non-) convergence of the perturbative series. Indeed, in principle, any converged physical quantity 
Q should be independent of the choice of values of these parameters. However, in practice, at a given 
loop order (L), all calculated physical quantities depend (artificially) on them: Q — >■ Q^^\a,b,a). 
Fixing a at the value obtained from the large order behavior, we consider that the optimal result 



43703-5 



B. Delamotte ef al. 



for Q at order L corresponds to the values (b)^. 
and a, that is, for which it is stationary: 



(L) (L) 



) of (6, a) for which Q depends the least on b 



Q 



(L) 
opt 



opt } 



with 



dQ^^^b, a) 
db 



dQ(^Hb,a) 
da 



^opt 



(2.18) 



where, of course, b^^l and a^j^l are functions of the order L. The validity of this procedure, known 
as the "Principle of Minimal Sensitivity" (PMS), requires that there is a unique pair (^optiQ^opt) 
such that Q^^^ is stationary. This is generically not the case: several stationary points are often 
found. A second principle allows us to "optimize" the results even in the case where there are several 
"optimal" values of b and a at a given order L: this is the so-called "Principle of Fastest Apparent 
Convergence" (PFAC). The idea underlying this principle is that when the numerical value of Q^^-* 
is almost converged (that is L is sufficiently large to achieve a prescribed accuracy) then the next 
order of approximation must consist only in a small change of this value: (5^^+^) ~ Q^^'. Thus, the 
preferred values of b and a should be the ones for which the difference between two successive orders 
Q,(i+i)) _ Q,(i)) is minimal. In practice, the two principles should be used 

together for consistency and, if there are several solutions to equation (|2.18|) at order L and/or 
L + 1, one should choose the couples {b^l,a'^l) and (&Qp^^\ a^p^^^^) for which the stationary 

values Q^^H^opt J ciopt) ^^^d Q^^'''^H^opt^^') Q^opt"^^) the closest, that is for which there is fastest 
apparent convergence. These principles have been developed and used in [23, 33, 31 , see also [lol |. 



3. 0(A/^)-symmetric theory 

Taking the above resummation procedure, let us show that it does not always eliminate un- 
physical FPs even in the simplest i/)^ model in d = 3. We study the /^u-function obtained within the 
massive scheme at four-, five- and six- loop orders [3l|. As was noticed in [27], the non-resummed 
/3„-function has no non-trivial root at the even orders of perturbation theory. Applying the re- 
summation procedure described by equations (|2.9l) - (|2.14p we find the Wilson-Fisher FP P close 
to the origin (see figure [5]) . There exists a large amount of studies of FP P and its properties are 
well-known (see e. g. ^^'^ references therein). However, we also find for larger w, in addition 

to P and for some values of a and 6, a new FP that we call S (in even orders of the loop expansion), 
see figure m Although this FP is unstable it changes the RG fiow structure and, if taken seriously, 
it would correspond to a tricritical FP (j2.1|l . However for the model under consideration the FP 
structure is well established [2] and this additional FP must be considered as an artifact of the 
fixed-d approach. 



y8„(u) 




Figure 2. Resummed /3„(u)-function of the 0(4) model at a = 6 and b — 4. The disks denote 
unstable FPs G and S, while the square corresponds to the stable FP P. 



Let us compare the stability properties of the critical exponents defined at P and S. Here, 
we consider u defined in equation (|2.4p . which is the leading correction to scaling exponent. We 
compute it at P and S from the resummed /3„-function and study its stability with the loop-order 
L and with respect to variations of the resummation parameters b and a. Note that although the 



43703-6 



Analysis of the 3d massive renormalization group perturbative expansions 



large order behavior of the series is known and leads to a preferred value of a j45|, [4^ it is possible 
to take it as a free parameter and to also study the stability of all the results w.r.t. variations of 
this parameter. We choose here to fix a, then to vary b and a so as to satisfy the PMS and PFAC 
and finally to check the stability of our results with respect to small variations of a. 

3.1. Wilson-Fisher FP 

First we consider uj for the Wilson-Fisher FP P. For 0{N) models, the analytically calculated 
value of a is known to be a = 0.14777422 x 9/(A^ + 8) ^] and we use it in our analysis. 




Figure 3. The exponent u) of the three-dimensional 0(4) model as a function of the resummation 
parameter h calculated at (a) FP P and (b) FP S. Results for FP P are presented at four- (grey 
dashed curve (green online)), five- (dark dashed curve (blue online)) and six-loop (solid curve) 
orders. The dots at each loop-order correspond to a stationary value of tj = uj{a,b) in both ct 
and b directions at the same time (L = 4: a = 3.2, 6 = 8.5, L = 5: a = 5.2, 6 = 9.5, L = 
a = 5, 6 = 13), at these values fastest apparent convergence between two successive loop orders 
is observed. Results for the FP S are presented in six and four loops at a = 6. No stationary 
point is found. The exponent u) is negative since S is repulsive. 

The results of our calculations are shown in figure [3] (a) where the exponent w for TV = 4 
(a = 0.110831) is shown as a function of the parameter h for the values of a for which stationarity 
is found for both b and a. The values obtained for four-, five- and six-loop orders are very close. 
An indicator of the quality of the convergence is given by the difference between the fifth and the 
sixth orders: ljj{L = 6) — ljj{L = 5) ~ 2.10^^. Our six loop estimation is a; ~ 0.783. This result is 
compatible with the results obtained by Guida and Zinn- Justin Q for = 4: a; = 0.774 ± 0.020 
(fixed d = 3, six/seven loops), w = 0.795 ± 0.030 (e-expansion, five loops), which is a check of the 
reliability of both the resummation scheme and the convergence criteria used here. 

3.2. Spurious FP 

We now analyze lo at the FP S. As expected, the results computed at S are very unstable with 
respect to variations of the parameters a and b. No stationary point is found for both b and a. To 
show the dependence of w on 6 at different loop orders we thus choose a given typical value of a, 
see figure [3] (b). A monotonous decrease of uj while increasing b is found. At fixed values of a and 
b we moreover find that increases with the loop order. 

At the end of this section some conclusion from the fixed d — ?> analysis of the 0{N) model can 
be made. Our results for w at P and S demonstrate two different features. For P, the convergence 
principles that we use — PMS and PFAC — allow us to find a reliable value of w, which is consistent 
with other estimations. On the contrary, w at S does not show any stationarity when b and a are 
varied. Furthermore, large differences between values in different loop orders (comparing to the 
results for FP P) are observed. 
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4. Cubic anisotropy 

The previous section was devoted to the study of the (j)'^ model with one coupling. We now 
consider a more complicated model with two couplings, namely describing cubic anisotropy. The 
effective Hamiltonian for this model reads Idl: 




The Hamiltonian (j4.1D is used to study the critical behavior of numerous magnetic and ferro- 
electric systems with appropriate order parameter symmetry (see e.g. [47|). The /3- functions are 
known at six-loop order in the massive scheme [48|. Their FPs structure is known The FP 
G (w* — V* — 0) and the Ising FP {u* — 0,v* ^ 0) are both unstable for all values of A''. Two 
other FPs: the 0{N) symmetric FP P {u* ^ 0,v* = 0) and the mixed one M {u* 0,v* ^ 0) 
interchange their stability at a critical value A^c of iV: for N < Nc the FP P is stable and M is 
unstable and vice versa for N > Nc- The critical value A^c has been found to be slightly less than 
3: for instance = 2.89 ± 0.04 in 01 and = 2.862 ± 0.005 in [13 . 




1.5 2.0 2.5 3.0 3.5 4.0 1.0 1.5 2.0 2.5 3.0 3.5 



u u 

(a) (b) 

Figure 4. Lines of zeroes of /?„ (solid curves) and (dashed curves) for the cubic model for (a) 
N = 2 with Q = 9 and 6 = 9 and (b) N = 3 for a = 8 and b = 10. The intersections of the solid 
and dashed lines correspond to different FPs: Wilson-Fisher P, cubic M and spurious Scub- The 
discs correspond to unstable FPs, while squares denote stable FPs. 

We have performed a stability analysis for the iV = 2, = 3 cases analogous to the 0{N) 
case. We take the 3d six-loop massive functions fS] and apply resumniation procedure generalized 
to the case of two variables (|2.15p - (l2.17p . Now, the FPs coordinates are defined by the system of 
equations: 

I3u{u,v) = 0, 

(4.2) 

l3y{u,v) = 0. 

Here, /3„ and /?„ stand for the resummed /3-functions of the model defined by (j4.ip . Being functions 
of the two variables u and v, Pu and /?„ describe surfaces in the space (u, v, j3). The FPs correspond 
to the points of common intersections of these two surfaces with the (m, u)-plane. As a guide for 
the eyes, we plot the lines of zeroes of the resummed jS^ and /?„ functions. The FPs, if they exist, 
correspond to the crossing points of such lines. To obtain such curves, one must use definite values 
of the resummation parameters a, &, a. As it has been shown in the parameter a for the cubic 



43703-8 



Analysis of the 3d massive renormalization group perturbative expansions 



model (|4J|) depends on the ratio z = v/u and is: a = 0.14777422 x (9/(2iV + 8) + z) for z > and 
a = 0.14777422 x {%/{2N + 8) + z/N) for -jff^ x < z < 0. The resummed /3-functions 

with the value a above are depicted in figure for N — 2 and iV = 3. 

From figure |4] one observes that, in addition to the usual FPs P and M, there exists another 
stable FP that has no counterpart within the e-expansion for both values of N, denoted by us as 
Scub- The presence of this FP, if taken seriously, would have important physical consequences since 
it would correspond to a second order phase transition with a new universality class. However no 
such transition has ever been reported. On the contrary, a first order behavior for all values of 
N larger than iVc and < is found within perturbative [3 S3 '^^ non-perturbative (s^ field 
theoretical analysis as well as numerical simulations [5l| in related systems (four-state antiferro- 
magnetic Potts model). Below we consider the dependence of physical quantities for the FP M 
and the spurious FP Scub upon variations of the resummation parameters. We then compare the 
results obtained with those for the 0{N) model. 

4.1. Mixed cubic FP 

Since the cubic model involves two coupling constants u and v the stability of its FPs is defined 
by two exponents, lji and U!2, which are the generalizations of equation (|2.4p for the case of two 
variables. At a stable FP, both lui and W2 have positive real parts. Calculating these exponents 




b 



Figure 5. The exponent uii at the mixed FP M for = 2 (a) and N = 3 (h) as a function of b 
at four (grey dashed curves (green online)), five (dark dashed curves (blue online)) and six (solid 
curves) loops. The dots at each loop order correspond to a stationary value of cj = ijj{a,b) in 
both a and b directions at the same time. For N — 2: a = 3.4, b — 14.3 {L — 4), a = 5.7, b — 8.3 
(L = 5), a = 5.8, b = 13.4 (L = 6); for TV = 3: a = 3.3, b = 14.6 (L = A), a = 5.3, b = 8.5 
{L = 5), a — 5.2, b — 12.6 (L = 6). At these values fastest apparent convergence between two 
successive loop orders is observed. 

for FP M (varying a and b in different L) we observe a similar picture as the one obtained for oj 
computed at P in 0{N) models (which is intuitively expected). In figure [5] we only present the 
behavior of the larger exponent loi as function of the resummation parameter. Stationary points in 
a and b are found. The values calculated at different loop orders at these points are very close to 
each other. For instance, for = 3 cj(L = 6) — uj{L = 5) w 0.001. Our six-loop results wi ~ 0.781, 
UJ2 = 0.011 for iV = 3 agree with other six-loop estimates uji = 0.781 ± 0.004, uj2 = 0.010 ± 0.004 
Q and LUi = 0.777 ±0.009, lu2 = 0.015 ±0.002 ^\E^. 

4.2. Spurious cubic FP 

We now study the FP Scub- First, we note that Scub for A = 2 and = 3 is an attractive 
focus, that is wi and uj2 are complex conjugate and Re(a;i)=Re(a;2) > 0. Our results are shown 
in figure |51 The real part of uii does not show any stationary values as function of resummation 
parameters (see figure [5]), and its behavior is similar to what was obtained for w at S in the 0{N) 
model. Moreover, the difference of values obtained in successive loop orders is always more than 
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0.3, that is 300 times more than for the FP M. It is thus clear that Scub is a spurious FP as is the 
FP S for the 0{N) modeL 



Re(6Ui) Re(6Ui) 




2 4 6 8 10 5 10 15 20 



(a) (b) 

Figure 6. The real part of exponent ui in the frustrated FP for A'^ = 2 at a = 5 (a) and A'' = 3 
at a = 6 (b) as a function of b at four (grey dashed curve (green online)), five (dark dashed 
curves (blue online)) and six (solid curves) loops. 



5. Frustrated magnets 

As in the previous section, we consider a model with two couplings but the one that describes 
frustrated antiferromagnets with non-collinear ordering. It is, for instance, used to describe a wide 
class of magnetic systems, such as antiferromagnets on stacked triangular lattice or helimagnets 
(see e.g. Contrarily to the previously described cases the physics of the phase transition in 

this model is not yet settled. For N = 2 and 3, different approaches predict a first order transition 
Q whereas it is found to be of second order in [T2|- The study of this problem within minimal 
subtraction scheme [55, already shed light on the reason of the discrepancy between these 
approaches showing that the stable FP should be considered spurious. We present here, within the 
massive RG scheme, additional arguments in this direction. 

The effective Hamiltonian relevant for frustrated systems is given by: 




where the (pi, i — 1,2, are A^-component vector fields and uq and vq are bare couplings that satisfy 
uq > and vq < 4uo ^ which corresponds to the region where the Hamiltonian is bounded from 
below. For > the ground state of Hamiltonian (|5.ip is given hy cpi = cj)2 = while for < 
it is given by a configuration where 4>i and 02 are orthogonal with the same norm. The Hamiltonian 
(|5.ip thus describes a symmetry breaking scheme between a disordered and an ordered phase where 
the 0{N) rotation group is broken down to 0{N — 2) (for details see |9| for instance). 

Let us first recall the picture of FPs obtained at leading order in e [M , [ssl - fssi . The FPs have 
two coordinates, u and v. For N larger than a critical value Nc{d), four FPs exist: apart from the 
usual Gaussian (u* — v* — 0) and 0{2N) [u* ^ 0,f* = 0) FPs, one finds the chiral-antichiral 
pair of FPs with coordinates u* > and v* > 0. The chiral FP C"*" is stable, whereas the other 
one, antichiral FP C~, is an unstable one. Above Nc{d), the transition is thus predicted to be of 
the second order for systems whose bare couplings lie in the basin of attraction of C+. When A^ 
is lowered at fixed d, C"*" and become closer and closer and finally collide and disappear at 
A^ = Nc{d). Below Nc{d), there is no longer any stable FP and the transition is expected to be of 
first order. The value of Nc{d) for d = 3 was a subject of intensive studies. In particular, it has 
been estimated within the e-expansion [13, , within pseudo e-expansion [ll| and within a NPRG 
approach [9|. A six-loop estimate has found Ac = 6.23 ±0.21 [Tlj within the pseudo-e expansion in 
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good agreement with the value Nc = 6.4 ±0.4 [501 derived from the resummed six-loops /3-functions 
computed in d — 3. According to this value of Nc, the physical systems that correspond to = 2 
or iV = 3 cannot undergo a second order phase transition. However a stable FP has been found for 
these values of TV in the fixed = 3 approach fl3\. Since this FP is not found in the e-expansion 
it is natural to wonder whether it is not an artifact of the fixed d approach as it is the case for 
the FPs S and S^uh in the 0{N) and cubic model, respectively. Thus, below we reconsider the 
/3- functions of the frustrated model (|5.1I) found in d = 3 at six loops |12| . Our analysis concludes 
that the stable FP found for frustrated systems at = 2, 3 i23| is spurious. 




0.0 0.5 1.0 1..5 2.0 0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 

u u u 

(a) (b) (c) 



Figure 7. Lines of zeros of Pu (solid curves) and fiv (dashed curves) for = 8 (a), = 6 (b) and 
A'^ = 4 (c) (30| . Results are obtained at a = 5 and b = 15. The crossings of solid and dashed lines 
correspond to different FPs: Gaussian (G), Wilson-Fisher (P), chiral (C^), antichiral (C~) and 
spurious (S"'", S~). 



First we investigate general situation with the change of A^. Parameter a for frustrated model 
has values depending on A^ [12]: a = 0.14777422 x 9/(2A^ + 8) for 4 (^jv^+s) > ^ > and a = 
0.14777422 x {9/{2N + 8) - z/2) for z < 0, z> 4,j2W+a) - Using this value in calculation we show 
the typical behavior (i.e. the behavior observed for the majority of a, b) of the lines of zeros of 
resummed /3-functions with fixed d = 3 in figure [7] for different A^, ranging from the large values 
of A^, A^ > A'c (figure [7] (a) to the small ones, N < (figure [7] (c). As one can see from the 
figures, the lines of zeros of /3„ form two branches, the upper and the lower one. Note that the 
lower branch has a form of a closed contour. This contour corresponds to zeros of function /?„ 
which exist already in the leading order in e. For N ^ 7 a. curve of zeros of /?„ intersects the closed 
contour generating the chiral-antichiral pair of FPs, C+ and C~ (see figure [7] (a)). When values 
of A" are between 7 and 6 this pair disappears in agreement with the scenario obtained in the 
e-expansion and with an estimate A^c ~ 6.23 [ll|. Then, till A^ « 5 there is no intersection between 
curves of zeros /3„ and (3^ in the region u > and v > and therefore there are no FPs in the 
given region (see figure [7] (b)). While at A^ < 5 we observe, that curve of zeros of I3y intersects the 
upper branch of zeros of /?„ forming two FPs, stable S"*" and unstable ones S~ (see figure [7] (c)). 
A similar situation has already been observed in [H^ that has led to a conclusion concerning the 
existence of another critical value of A^. From the above sketched analysis one can see that the FP 
solutions C"*", C~, on the one hand and S+, S~ on the other hand correspond to the crossing of 
different branches of the /?„ curve. In the first case (C"*", C~) this is the lower branch, whereas in 
the second case (S^, S~) this is the upper one. Below, we will show the similarities in the behavior 
of solutions S+, S~ to those that were obtained in the former section |3] for the spurious FP. In 
the following subsections we separate the analysis of stable FPs for cases "large A^" and "small A^" 
meaning by this two different regimes when A^ is greater than A"c (c.f. figure[7] (a)) or smaller than 
A'c (figure [7] (c)), correspondingly. 
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5.1. Case of large TV 

Similarly to the above analyzed cubic model here we present the results of the analysis for 
stability exponents at d = 3. We start with the data obtained for the stable FP C+ at iV > iVc. 
In this case we take N = 8, therefore a = 0.0554, and analyze how the exponent w changes with 
a and b. We find that the PMS is satisfied at four-, five- and six-loop orders: for suitable values of 
the parameters b and a the two exponents uji and lu2 depend weakly on these parameters and are 
reasonably well converged. This is clear from figure |51 where we show the 6-dependence of wi and 
W2- Moreover, the difference between the values at five and six loops of, for instance, lj2 is small: 
UJ2{L = 6) —uj2{L = 5) ~ 0.012. Note that our values of wi and cj2 in this case are fully compatible 
with those obtained in |56] . 

0)2 
0.39 F 




0.820 
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0.37 \ 
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0.35 ; 

0.34: 

0.33 : 
0.32^ 
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(b) 



Figure 8. The exponents aji(a) and uu2 (b) at FP of frustrated model for JV = 8 as a function 
of b at four (grey dashed curves (green online)), five (dark dashed curves (blue online)) and six 
(solid curves) loops. The dots at each loop-order correspond to a stationary value of uj — cj(a, b) 
in both Q and b directions at the same time. For ui: a = 5.1, b = 14.8 (L = 4), a = 7.8, b = 17 
(L = 5), a = 7, 6 = 12.2 (L = 6); for uj2: a = 6.5, b = 27.1 (L = 4), a = 7.9, 6 = 16.4 (L = 5), 
a = 7.5, b = 11.7 (L = 6). 



These results indicate that the convergence properties of the iV = 8 frustrated model are 
globally similar to those of the physical FPs of the 0{N) and cubic models although very likely 
less accurate because in the latter case the resummation is less efhcient due to the presence of 
complex symmetry. 



5.2. Case of small N 

Now, let us study the peculiarities of the stable FP S+ in the region u > 0, v > for the 
physical cases A*" = 2, 3. The results obtained at the calculation of w at this FP with variations for 
b and a are given in figure [5] 

As it was noted in stability exponents attain (in most cases) complex values at FP S"*" with 
positive real parts of uj indica ting a stable focus topology. We, again, take for a the value obtained 
from the large order analysis [13|: a = 0.1108 for = 2 and a = 0.095 for iV = 3. For these values 
of a and for L = 4,5 and 6, we find that Re(wi) (or equivalently Re(w2)) considered as a function 
of a and b is nowhere stationary, even approximately, see figure Moreover, at fixed a and b, 
the gap between the values of Re((u;i) at two successive loop-orders: Ke{uji){L -1-1)— Re(a;i)(L), is 
always large, of order 0.5 for iV = 2 and 0.2 for = 3, see figure |H1 Thus neither the PMS nor the 
PFAC are satisfied for these values of N. 

The 6-dependence of uji obtained in figure shows the similarity with the results for the FP S 
of the 3d 0{N) model (compare with figure [3] (b)) as well as with Scub of cubic model (compare 
with figure in]). From the figure [2] the bad convergence properties for spurious FP are evident. Also 
the value of oj in FP S+ decreases when increasing b, just as it was observed for the FP S and 

Scub- 



43703-12 



Analysis of the 3d massive renormalization group perturbative expansions 



Re(wi) Re(ctJi) 




b 



Figure 9. The real part of the exponent coi in the FP S"*" for = 2 (a) and = 3 (b) as a 
function of b for a = 6 at four (grey dashed curves (green online)), five (dark dashed curves 
(blue online)) and six (solid curves) loops. 



Thus we observe a situation analogous to the previously studied models. Similarly to the phys- 
ical FPs of these models, for FP C+ at > A'c, we observe good convergent behavior, which is 
displayed via weak change of exponents and small difference between the results of different loop 
order. While for FP S+ at < A^c we do not observe the regions of stabilization of critical expo- 
nents. The values of stability exponents considerably differ in different orders. Since such properties 
are characteristic of FP S and S^uh, this supports the unphysical character of FP S+. 



6. Conclusions 

The main motivation of the study reported in this paper was to attract attention to some 
problems that arise in the course of RG analysis of critical behavior of complex Hamiltonians. 
Being formulated briefly one can state that RG has proven its accuracy when it is used to calculate 
numerical values of different universal observables that govern the second order phase transition and 
critical behavior in more general sense. However, sometimes RG is used to prove the occurence of a 
critical point itself, fn such cases, the occurence of stable and reachable FP of a RG transformation 
is considered to be the proof of a second order phase transition. This is the problem we addressed in 
our paper. Namely, given that in principle, a non-physical solution may arise in any computational 
scheme, how could one single-out such solutions from the physical ones? 

In particular, several different solutions for FP may be obtained in the analysis of the expressions 
of /3-functions obtained within the perturbative RG as series in coupling constant without use of 
e- or pseudo-e expansions (27l . [s^ directly at fixed space dimension d. It is clear that the number 
of such solutions increases with an increase of the loop approximation. The way to distinguish 
between physical and unphysical FPs is resummation procedure, which is expected to eliminate all 
spurious solutions. However, generally this is not the case. As it follows from the analysis performed 
in this paper for the 0{N) model, except physical Wilson-Fisher FP, another FP survives. The 
current knowledge of the properties of 0{N) model at d = 3 excludes the occurence of such 
FP. Additionally, the obtained FPs clearly differ in their convergent properties. In particular, the 
behavior of spurious FP is strongly dependent on the fit parameters used in resummation and does 
not satisfy the principles of convergence (PMS, PFAC). A completely identical picture is observed 
for the cubic model, whose critical behavior is also known. New (spurious) FP, found for this model 
by a fixed d approach, demonstrates the same behavior as the spurious FP of 0{N) model. 

Although the problem of discerning between physical and unphysical FPs is quite general. Here 
we address it taking as an example the phase transition for the model of frustrated magnets. The 
question of order of phase transition in such systems remains unclear. The experimental results 
make it possible to divide the conditionally tested materials into two groups that differ by their 
critical exponents. The feature of one group is negative exponent 77, while the scaling relations are 
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violated for the critical exponents of the other group (see review in [^). Both phenomena should 
not occur at the second order phase transition. Recent numerical simulations give evidence of the 
first order phase transition for frustrated systems (57l - l6l| . The RG investigations based on the 
e- and pseudo-e-expansions [ll, 1^ 17 1 indicate the absence of the second order phase transition 
within (j>* theory for a frustrated model at physical values N = 2 and N = 3. This was corroborated 
within NPRG Sm-ill. However, in the RG analysis at fixed d with application of resummation 
procedures a stable FP was found at iV = 2, 3, but only starting from four-loop approximation 
|l2fl . Since the topology of this FP is stable focus, the difference in the experimental data was 
explained by an approach of RG flows to this focus. 

The properties of such FP were already checked within the minimal subtraction scheme [2^ , 
which makes it possible to continuously follow the behavior of the FPs with the change of space 
dimension from d = 3 to d = 4, where only Gaussian FP describes the physics of a system. There, 
it was shown that the spurious FP found at d = 3 persists at d = 4 as well. Such an observation 
enabled us to assume that this FP is an unphysical one. 

In the present study the expressions obtained in the massive scheme at d = 3 are used. There- 
fore, it is impossible to trace them to d = 4 by continuous change of d. However, the analysis we 
performed here shows that stable FPs obtained for values N = 2, 3 on the one hand and for the 
values higher than the critical value one the other hand, differ in their behavior in respect to 
a variation of fit parameters of the resummation procedure. While for values oi N > Nc reliable 
results may be chosen by principles of convergence, for N < el strong dependence of FPs be- 
havior on the resummation parameters does not make it possible to apply these principles. The 
latter fact serves as an argument to consider FPs for N < to be unphysical ones. 
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AHa/ii3 3d Maci/iBHi/ix posBUHeHb nepTypGarnBHoY 
peHopMa/ii3ai4iMHoV rpyni/i: fle/iiKarHa cnpasa 

B. fle;i5iMo-iEl, M. flyflK^, K). To/ioBaJ^, fl,. MyaH^l 

JlaGopaxopia xeopexuHHoT 4>i3MKi/i KOHfleHCOBaHi/ix cepeflOBi/im, CNRS-UMR 7600, VHiBepcuTeT fl'epa i 
MapiV Kropi, 75252 napM>K Cedex 05, 0paHL4iiR 

iHCTMTyx 4>i3MKM KOHfleHCOBa HMX CMCTeM HAH VKpaYHi/i, UA— 79011 JlbBie, VKpaiHa 
iHCTMxyT TeopeTMHHoT cfji^MKM, VHiBepci/iTeT MoraHa Ken/iepa, A— 4040 /IImli, AscTpia 

AHajiisyeTbca ecj^eKTHBHicTb nepTyp6aTHBHoro niflxofly y Mejofli peHopMajiisaLiiMHoT rpyni/i npi/i 4>iKCOBaHiH 
BMMipHocxi npocTopy d B xeopiT kpmtmhhmx aBwm. Posmaflaexbcn jpi/i MOfle.ni: 0{N), Ky6iHHa i aHTi/icfje- 
poMarHixHa MOflejib OBHaneHa Ha Tpi/iKyxHiM rpaTL4i. Mi/i po3maflaeMO yci MOflejii npw ctiiKCOBaHOMy d = 3 \ 
aHa/iisyeMO npoL4eflypy nepecyMOByBaHHa, aKa Bi/iKopi/icTOByeTbca fljia o6Hnc/ieHHa kphtmhhmx noKaBHMKis. 
Ml/1 cnoHaTKy noKasyeMO, mo fl/ia 0{N) MOflejii nepecyMosyBaHHa He BMKJiKDHae yci necJjiai/iHHi (4>iKTHBHi) 
HepyxoMi TOHKM (HT). Tofli yBa>KHO BUBHaexbca 3a/ie)KHicTb cfjiKTHBHoT HT, a TaKO>K HT Bi/ibcoHa-Oiiuepa 
Bifl napaMexpiB nepecyMOByBann. Kpi/iTi/iHHi noKasHi/iKi/i b HT Bi/ibcoHa-Oiiuepa c/ia6KO 3ajie>KaTb sifl nape- 
MexpiB nepecyMOByBaHHa. Ha npoxi/iBary L4bOMy, noKaBHi/iKi/i b cfjiKTMBHiM HT, a TaKO>K caMe TT icnyBaHHa, 
cyxTeBO 3ajie>KaTb bifl limx napaMexpiB. KySiMHoT MOflejii oxpuMaHO Hosy cxiwKy HT i noKasaHO, lu,o i TT 
B/iacxuBocxi xaKO>K BHaHHO 3ajie>Kaxb sifl napaMexpiB nepecyMOBysaHHa. BoHa BwaBJiaexbca ctiiKxi/iBHOK), jik 
oniKyBajiocb. U-loflo cfjpycxpoBaHoT MOflejii, xo icHyraxb flBa Bi/inaflKi/i b 3a/ie>KHOcxi sifl BHaneHHa Hucjia cni- 
HOBMX KOMnoHeHX. KojiM 6ijibLue 3a KpuxMHHe BHaneHHiR A^c, noBefliHKa y cxiMKiw HepyxoMiCi xoHL4i nofli6Ha 
flo noBefliHKM y HT Bi/ibcoHa-0iiuepa. Ha npoxuBary L4bOMy, fl^na A' < A^c pesyjibxaxn oxpuMaHi fljia cxitiKoT 
HT nofli6Hi flo xwx, mo oxpuMaHi b 4>iKXWBHi/ix HT 0{N) i Ky6iHHoT MOflejii. 3 L4boro anamsy mm po5hmo 
BMCHOBOK, LU,o cxiwKa HT BHawfleHa fljia A' < A^c b cfspycxpoBaHiw MOfle/ii, e 4>iKXMBHOio, OcKijibKi/i A^c > 3, 
MM po6mmo BMCHOBOK, LU,o nepexifl fljiji XY i raH3eH6epriBCbKoro cjjpycxpoBaHoro MarnexHKa e cfiaBOBMM 
nepexoflOM nepiuoro pofly. 

K/iiOHOBi c/iOBa: Teopia no/ifi, peHopMam3au,it^Ha rpyna, Kpi/iTi/iHHi aBi/ima, reopifi 36ypeHb, 
nepecyMOByBaHHfi 
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